Self-energy-functional approach to systems of correlated electrons 
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The grand potential of a system of interacting electrons is considered as a stationary point of a 
self-energy functional. It is shown that a rigorous evaluation of the functional is possible for self- 
energies that are representable within a certain reference system. The variational scheme allows 
to construct new non-perturbative and thermodynamically consistent approximations. Numerical 
results illustrate the practicability of the method. 
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I. INTRODUCTION 

Systems of strongly correlated electrons continue to 
represent a central subject of current research. Dif- 
ferent interesting correlation phenomena, such as high- 
temperature superconductivity^ Mott metal-insulator 
transitions^ or itinerant ferromagnetism^ are far from 
being finally clarified. Progress in this field crucially de- 
pends on the development of new theoretical methods as 
even highly idealized model systems pose notoriously dif- 
ficult problems. There arc only a few general approaches 
which are able to access the equilibrium thermodynamics 
as well as excitation properties of an extended system of 
correlated electrons. 

General methods can be based on the Green's-function 
formalism of Luttinger and Ward^ and Baym and 
Kadanoff£ Here the grand potential fl is expressed in 
terms of the time- or frequency-dependent one-electron 
Green's function G. The functional 0[G] can be shown 
to be stationary at the physical G. In principle, this is 
an exact variational approach which provides informa- 
tion not only on static equilibrium but also on dynamic 
excitation properties. The functional dependence S1[G], 
however, is generally not known explicitly as it must be 
constructed by summation of an infinite series of renor- 
malized skeleton diagrams. In the standard approxima- 
tion the exact but unknown functional is replaced by an 
explicitly known but approximate one which is based on 
an incomplete summation of the diagram series. This 
leads to the well-known perturbational ("conserving") 
theories^ Higher-order approaches as the fluctuation- 
exchange approximation^ are mainly applied to discrete 
lattice models while for continuum systems, e.g. for the 
inhomogeneous electron gas, one has to be content with 
lowest-order theories as the GW methodi 7 ! 8 ' 9 

A second type of general methods is based on density- 
functional (DF) approaches piSiii Normally these aim at 
the inhomogeneous electron gas but can also be applied 
to Hubbard-type lattice models^ Compared with the 
Green's-function formalism, there is a conceptually simi- 
lar situation for DF approaches: In the latter the ground- 
state energy E (or the grand potential is given as 
a functional of the (static) density n. The variational 
principle associated with the functional E[n] is rigorous 
but cannot be evaluated as E[n] is generally unknown. 



In the standard local-density approximation (LDA) the 
(unknown part of the) functional is replaced by an ex- 
plicitly known but approximate functional which is taken 
from the homogeneous system. For systems with weakly 
varying density the LDA should be justified. Informa- 
tion on excitation properties is contained in dynamic 
response functions which are in principle accessible via 
time-dependent DF theory^ where the action A is con- 
sidered as a functional of the time-dependent density n. 
Again, the exact but unknown functional A[n] is approx- 
imated to make it explicit and the variational principle 
is exploited afterwards. 

The method proposed here rests on a variational prin- 
ciple which uses the electron self-energy £ as the basic 
dynamic variable. A new functional f2[S] is constructed 
which can be shown to be stationary at the physical self- 
energy. The main result is that the variational principle 
can be exploited without any approximation of the func- 
tional dependence. Namely, a rigorous evaluation of the 
functional f2[S] is possible on a certain subspace of trial 
self-energies. Trial self-energies must be representable 
within an exactly solvable reference system sharing the 
same interaction with the original system. 

This result has important consequences as it opens 
a route for constructing a novel class of approxima- 
tions. Although the self-energy essentially contains the 
same information as the Green's function or the (time- 
dependent) density, the new approach is conceptually 
contrary to the Green's-function approach and to the 
DF approach as there is no approximation to be tol- 
erated for the central functional. Instead of approx- 
imating the functional itself, it is considered on a re- 
stricted domain. The self-energy-functional approach is 
completely general and yields approximations which are 
non-perturbative, thermodynamically consistent and sys- 
tematic. Opposed to numerical techniques directly ap- 
plied to systems of finite size, the self-energy-functional 
approach provides a variational or self-consistent embed- 
ding of finite systems and thus yields results in the ther- 
modynamical limit. Such techniques are needed to con- 
struct phase diagrams from standard correlated lattice 
models. A potentially fruitful field of application are sys- 
tems with competing types of order resulting from spin, 
charge or orbital correlations as it is typical e.g. for nu- 
merous transition-metal oxidesi^ 
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In the present paper the approach is introduced and 
a number of general aspects are discussed in detail (Sec. 
ITTj) . To demonstrate its usefulness, two applications will 
be considered for the single-band Hubbard model: In 
Sec. IIIII it is shown that the dynamical mean-field the- 
ory (DMFT)^ can be recovered within the self-energy- 
functional approach, namely by choosing a decoupled set 
of impurity Anderson models as a reference system. The 
DMFT generally requires the treatment of a quantum- 
impurity problem with an infinite number of degrees of 
freedom (n s = oo). In Sec. II VI a new approximation is 
discussed which is based on an impurity model with a 
finite number of degrees of freedom only and which ap- 
proaches the DMFT for n s = oo. The method is closely 
related to the cxact-diagonalization approach (ED)^£»i£ 
Opposed to the ED, however, thermodynamical consis- 
tency is guaranteed at any stage of the approximation. 
New approaches beyond the mean-field level will be dis- 
cussed elsewhere. The conclusions and an outlook are 
given in Sec. 

II. SELF-ENERGY-FUNCTIONAL APPROACH 

Consider a general Hamiltonian H = Ho(t) + i?i(U) 
with one-particle ("hopping") parameters t and two- 
particle interaction parameters U: 

H = Y1 t ^a c fi + o U a ^clc^CjC S . (1) 

Here a, /3, ... refer to an orthonormal and complete set of 
one-particle basis states. We are interested in the equi- 
librium thermodynamics and in elementary one-particle 
excitations of the system for temperature T and chem- 
ical potential (i. This is described by the one-particle 
Green's function G a p{vJ) = ((c Q ;c^}} of the imaginary 
fermionic Matsubara frequencies iuo = i(2n + 1)ttT with 
integer The Green's function can be calculated from 
the self-energy T, a p(iuj) via the Dyson equation. Us- 
ing a matrix notation, this reads as G = Go + GoSG 
where Go = l/{iu> + [i — t) is the "free" Green's function. 
The self-energy is given by S = S[G] = T^S^Gl/SG, 
where $[G] is the so-called Luttingcr-Ward functional!^ 
This allows to derive the Green's function from a varia- 
tional principle: One has <5fi[G]/<5G = where fi[G] = 
$[G] + Trhi(-G) - Tr((Gg 1 - G^G) and using the 
notation TrA = A aa (iuj). In general, however, 

the functional $[G] is not known explicitly which pre- 
vents an evaluation of fi[G] for a given G. So-called 
conserving approximations^ provide an explicit but ap- 
proximate functional $ pcrt .[G] ss $[G]. However, these 
are weak-coupling approaches where a certain subclass of 
$ diagrams is summed up. 

Here a different but still rigorous variational principle 
is proposed which is based on a functional G = G[53] 
defined as the inverse of 53 = 53 [G]. We can assume the 
latter to be invertible (locally) provided that the system 



is not at a critical point for a phase transition (see Ap- 
pendix ^J. Consider then: 

fi t [53] = TrlnHGo 1 - 53)" 1 ) + F[E] (2) 

where F[E] = $[G[53]] - Tr(53G[53]) is the Legendre 
transform of 3>[G]. The subscript t indicates the explicit 
t dependence of fi due to the free Green's function Go- 
Using T^SFm/SIl = G[53], one finds that 

<Jft t [E]/<SE = <s> G[53] = (G- 1 - E)- 1 . (3) 

Thus fit [53] is stationary at the exact (physical) self- 
energy and its value is the exact grand potential of the 
system. Again, the problem is that the functional fit [53] 
is in general not known explicitly. 

As the domain of the self-energy functional fit [53] we 
define the class of all t' rcprcscntable self-energies. 53 
is termed t' representable, if there is a set of hopping 
parameters t' such that 53 is the exact self-energy of 
the model Ho(t') + Hi(U). This implies that any self- 
energy in the domain of fit [53] can be parameterized as 
53 = 53 (t'). The interaction parameters U are taken to 
be fixed. Suppose we are interested in the model H = 
H (t) + ffi(U). Then the function fi t (t') = fit [53(f)] is 
stationary at t' = t. Thus dQ t (t')/dt' = 0. 

It is important to note that F[S] is universal: The 
functional dependence is the same for any t, i.e. it re- 
mains unchanged for an arbitrary reference system H' 
with the same interaction but modified hopping param- 
eters: H' = H (t') + ffi(U). F[E] is universal as it is 
the Legendre tranform of 3>[G] which in turn is univer- 
sal because it can be constructed formally as the sum 
of all closed, irreducible, and renormalizcd skeleton dia- 
grams which, apart from G, include the vertices U only. 
Consequently, one has: 

Q t ,[E] = TrlnHGo- 1 - 53)- 1 ) + F[53] , (4) 

for the reference system H' with G 1 = iu) + fi — t'. 
Combining Eqs. and J3J, F[S\ can be eliminated: 

fi t [53] = nt,[E] + TrlnHG^-E)- 1 ) 

- TrlnHGo- 1 -E)- 1 ). (5) 

Evaluating the functional fit [53] for self-energies param- 
eterized as 53 = 53 (t'), one obtains: 

fit [53(f)] = f2' + Trln(-(G- 1 -S(t')) -1 ) -Trln(-G') . 

(6) 

Here it has been used that fit' [S(t')] = fi', the ex- 
act grand potential of the reference system H', and 
(G 1 — 53(f)) -1 = G', the exact Green's function of 
W . Suppose that the reference system H' is much sim- 
pler than the original system H so that it can be solved 
exactly for any t' belonging to a certain subspace of the 
entire space of hopping parameters. The resulting Eq. © 
is remarkable, as it shows that the functional fit [53] can 
be evaluated rigorously for trial self-energies 53 = 53(f) 
taken from the reference system H'. 
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This is the main result. Contrary to previous ap- 
proaches (e.g. conserving theories, LDA), there is no need 
to approximate the functional dependence in a funda- 
mental variational principle. Approximations are con- 
structed by searching for a stationary point of Qt [S] on 
a restricted set of trial self-energies S(t'). 

The stationary point is determined by the Euler equa- 
tion: dil t [E(t')]/dt' = 0. Calculating the derivative, 
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FIG. 1: Schematic representation of the Hubbard model H 
(a), an equivalent model (b), and a possible reference system 
H' (c). See text for discussion. 



Note that the equation involves, apart from Go, quanti- 
ties of the reference system H' only. The linear response 
of the self-energy of H' due to a change of the hopping t' 
can be calculated along the lines of Ref. 0. It turns out 
that d~E(t')/dt' is given by a two-particle Green's func- 
tion of H'. Since G' = G[S(t')], the exact self-energy 
of the system H is determined by the condition that the 
bracket in (7J| be zero. Hence, one can consider Eq. Q 
to be obtained from the exact equation that determines 
the "vector" S in the self-energy space through projec- 
tion onto the hypersurface of t' representable trial self- 
energies S(t') by taking the scalar product with vectors 
dH(t')/dt' tangential to the hypersurface. 

An analysis of the second derivative 
d 2 Qt[£(t')]/dt' a gdt' s shows that a stationary point 
is not an extremum point in general. This feature is 
shared with the time-dependent DF approach^ the 
Green's-function approach^ and also with a recently 
considered variant^ Only in the static DF theory there 
is a convex (density) functional>i2*ii*i£ Nevertheless, the 
proposed self-energy- functional approach is systematic: 
For any sequence of reference systems H' including 
more and more degrees of freedom and converging to 
the original system H there is, from the variational 
principle, a corresponding sequence of grand potentials 
which must converge to the exact ft = Sl t [£(t)] as the 
subspace of trial self-energies increases and eventually 
includes the exact self-energy S(t). 

III. RELATION TO THE DMFT 

Given an original model H, what could a suitable ref- 
erence system H' look like? Consider, for example, H 
to be the Hubbard models which is shown in Fig. la 
schematically: A hlled dot represents a correlated site 
i with on-site Hubbard interaction U, and a line con- 
necting two sites i and j represents the nearest-neighbor 
hopping tij. The number of sites is L i— > oo. Fig. lc 
shows a conceivable reference system H'. H' is obtained 
from H (Fig. la) by (i) adding to each correlated site i a 
number of n s — 1 uncorrelated ("bath") sites k = 2, n s 
(open dots) which are disconnected from the rest of the 
system, by (ii) switching off the hopping tij between the 
correlated sites and (hi) switching on a hopping V^jt to 
the bath sites. After step (i) the Hamiltonian Fig. lb 



(in the figure n s = 5) has an enlarged Hilbert space 
but the same self-energy. It is important to note that 
steps (i) - (hi) leave the interaction part unchanged and 
thus preserve the functional dependence F[E}. Actually, 
the system H' is a set of L decoupled single-impurity 
Anderson models (SIAM)2i with n s sites each. Com- 
pared to H, the problem posed by H' is strongly sim- 
plified. This is achieved at the cost of restricting the 
set of trial self-energies. In particular, as the correlated 
sites are decoupled in H', the trial self-energies are local: 
Hij(iuj,t') oc Sij. One has to consider H' for arbitrary 
one-particle parameters, namely the on-site energies of 
the correlated ("c") and of the bath sites ("a"), and 

efl, respectively, and the hopping ("hybridization") Vi t k 
between them and take these as variational parameters 
in the principle <5f2t[X(t')] = 0. 

Let us discuss the case n s i— > oo. For a homoge- 
neous phase of the (translationally invariant) original 
system, f2t[E(t')] will be stationary at a homogeneous 
set of variational parameters: t' = {e\ c \ e\ \ J, V^fc} = 

{e(°\ e;° , Vfe}. Consequently, it is sufficient to consider 
one SlAM only. As the different equivalent SIAM's are 
spatially decoupled, not only the self-energy but also its 
linear response is local: 9Sy (t') / dt' oc 5^. To solve the 
Euler equation Q), it is thus sufficient to fulfill the "lo- 
cally projected" equation 

( r-i r . r vr s ) =G' tl M- (8) 
\G (iw) - H{iw)J u 

This is just the self-consistency equation of the DMFT: 15 
The SIAM parameters have to be found such that the on- 
site ("impurity") Green's function at a correlated site i 
coincides with the on-site Green's function of the Hub- 
bard model which is calculated from Go and the ( "impu- 
rity") self-energy of H' by means of the Dyson equation. 
Therefore, one can state that the DMFT (as an approxi- 
mation for any finite-dimensional system or as the exact 
theory in infinite dimensions) is recovered as a station- 
ary point of f2t[E] when restricting the search to local 
self-energies representable by a SIAM. 

Within the DMFT the computation of the self-energy 
requires an iterative procedure: S id >— ► S ncw . Here it 
turns out that this corresponds to a certain (discrete) 
path on the hypersurface of SIAM trial self-energies. 
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Convergence of the iteration (S id = S ncw ), however, is 
by no means guaranteed physically but depends on the 
contracting properties of the map S id >— > S ncw . The 
self-energy-functional approach offers an alternative as 
instead of solving Eq. © one may calculate f2t[E(t')] 
by Eq. (j^J and determine the stationary point. Hence, 
the DMFT can also be obtained by a more direct compu- 
tation avoiding any iterations - similar (in this respect) 
to the random-dispersion approximation^ Note that in 
case of more than a single stationary point there is also 
an equally direct access to metastablc phases. 

For any inhomogeneous situation, Eq. (jHJ represents 
a system of self-consistency equations to fix the param- 
eters of non-equivalent impurity models labeled by the 
site index i. The models can be solved independently 
but are coupled indirectly due to the matrix inversion 
in JHJ. This exactly recovers the DMFT generalized to 
systems with reduced translational symmetryi2i24 
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FIG. 2: Grand potential Q (per lattice site) and the differ- 
ent contributions (per lattice site) according to Eq. @: Q' , 
Trln(-(Go 1 - S(t'))" 1 ), and Trln(-G') for U = 4, T = 0, 
and \i = U/2 (half-filling) as functions of V . Stationary points 
(arrows) at V = ±0.519 (metal) and at V = (insulator). 
e c — and e a = 2. 



IV. A CONSISTENT ED METHOD 

A brief discussion of two limiting cases of the Hubbard 
model may be instructive. Consider the band limit with 
U = first. Here H = J^ija ^ij c \<y c j<y describes a system 
of non-interacting electrons. This case is exceptional as 
obviously the functional F[E] = and therefore fit[X] = 
Trln(— (Gq 1 — E) _1 ). Any valid reference system H' 
must have the same (i.e. a vanishing) interaction part as 
H, and thus E(t') = and ft t [£(t')] = T r M- G o" 1 )> 
the exact grand potential for non- interacting electrons. 

The atomic limit, H = ^^{to^Cia + {U /2)n i(T m- a ) 
is more interesting as $[G] and F[S\ cannot be con- 
structed explicitly. Within the self-energy-functional ap- 
proach one has to compute S(t'), G', and CI' for a suit- 
able reference system H' and to insert into Eq. (JSJ) for 
optimization. The only meaningful choice for the ref- 
erence system is H' = H in this case. Obviously, this 
yields the exact solution. Generally, whenever the origi- 
nal model H is exactly solvable, the choice H' = H will 
do. 

For a non-trivial model H, the choice H' = H is use- 
less for a practical computation. Any simplified reference 
system, however, yields a consistent approximation. The 
case of the DMFT discussed in Sec. Illll is an illustrative 
example. On the other hand, in the context of the DMFT 
actually both, H and H' , are highly non-trivial models, 
and further approximations or large-scale numerics are 
needed to treat the reference system H' . More simple 
approximations for the Hubbard model which are still 
constistcnt are generated by considering reference sys- 
tems with a finite number of degrees of freedom. The 
reference system of Fig. lc with n s < oo is an interest- 
ing example which shall be discussed in the following. 
For small n s one can easily obtain numerical results as a 
complete diagonalization of H' is feasible. 



Consider the Hubbard model 

H = Yl *« c L c jV + n la n^ a (9) 

at temperature T = and chemical potential /i = U/2. 
For symmetric conditions this implies half-filling. The 
Hamiltonian of the reference system is given by H' = 

Ei^'W with 

H'(i) = ^ ^Aa^o + I} n ™ n i-° 
a a 

+ e< tk a \ka a ik<y + Y ( V ^k4 a a ika + h.c.J . 

(7,k= 2 a,k 

(10) 

For the sake of simplicity we consider a homogeneous 
paramagnetic phase and the most simple case n s = 2 
where one is left with three independent variational pa- 
rameters only, namely V = Vi.k=2, (-a = e^^? =2 , and 

e c = e[ c \ The computation of the different contributions 
to the grand potential, Eq. is straightforward: Diag- 
onalization of H' yields the ground-state energy Eq and 
O' = Eq — n(N') as well as the excitation energies, the 
ground state and the excited states. The Green's func- 
tion G' and the free Green's function Gg can be com- 
puted from their respective Lehmann representations. 
The self-energy of the reference system is obtained as 
S(t') = G 1 — G'~ . Since the self-energy is local, 
as in the DMFT, the lattice structure enters via the 
free (U = 0) density of states only. Therefore, the It- 
sum which appears in the first trace in Eq. © can be 
performed conveniently by a one-dimensional density-of- 
states integration. A semi-elliptical free density of states 
with the band width W = 4 is used for the calculations. 
This sets the energy scale for the results discussed below. 
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Fig. |21 shows the grand potential fl and the three dif- 
ferent contributions as functions of V. The interaction 
is kept fixed at U = W = 4. The remaining variational 
parameters are set to e c = and e a = U/2 = 2, as 
required by particle-hole symmetry. Each of the three 
contributions strongly depends on V and none of them 
has a stationary point at a finite V =/= 0. Two of them 
show a singular behavior at V — 0. Contrary, the result- 
ing f2 is regular for any V and shows a much weaker V 
dependence. There are three stationary points which are 
indicated by the arrows. The maximum at V = corre- 
sponds to an insulating phase since for n s = 2 and 
V = is the Hubbard-I self-energy which implies a van- 
ishing spectral density — (l/7r)ImG(w + i0 + ) at u = 0. 
The minima at V = ±0.519 correspond to a metallic 
phase. f2 as well as the different contributions are sym- 
metric functions of V. As S(VJw) = E(— V,ui), however, 
this symmetry is trivial and docs not yield an additional 
physical phase (see also Appendix^) . Due to the lower f2 
the metallic phase is stable as compared to the insulating 
one. 

The minimum at V = 0.519 is actually a saddle point 
if the entire space of variational parameters is consid- 
ered. This is demonstrated by Fig. which shows 
as a function of e c for fixed V = 0.519 and e a = 2. 
While fl(V) is at a minimum for V = 0.519, f2(e c ) is 
at a maximum for e c = 0. In the (V, e c ) space one 
therefore encounters a saddle point. As already noted 
in Sec. EI there is no reason to expect an extremum 
in general. It is worth mentioning that stationarity at 
e c = is consistent with the requirements of particle- 
hole symmetry. For any e c ^ the impurity model is 
asymmetric. This can be seen from the inset where the 
average occupations of the impurity and of the bath site 
are plotted as functions of e c . The total particle number 

( N ') = Ea«4<V> + (4<V)) = n c + n & (i and k = 2 
fixed) is constant: (N 1 ) = 2. 



With respect to the third variational parameter e Q , the 
grand potential SI is at a maximum for e a = 2 = fx, see 
Fig- El Again, this value is required by particle-hole sym- 
metry. If e a exceeds a certain critical value (away from 
the stationary point), the ground-state of the reference 
system H' no longer lies within N' — 2 subspace but 
is found in the N' = 1 or N' = 3 subspace, respec- 
tively (see inset of Fig.0J. While £1' is continuous at the 
level crossing, the symmetry of the ground state changes. 
Consequently, there is a discontinuous change of the trial 
self-energy which implies a discontinuous change of Q. 

Consider now the original model at slightly modified 
parameters, e.g. U, /x, or T. Clearly, the stationary point 
of il will be expected then at slightly different values of 
the variational parameters V, e c , e a . This implies that 
all physical quantities which derive from the thermody- 
namical potential f2 will be continuous functions of the 
(original) model parameters in general - irrespective of 
the fact that the reference system includes a finite num- 
ber of degrees of freedom only: It is a typical feature 
of any mean-field approach that results are directly pro- 
vided for the thermodynamical limit. A discontinuous 
jump of f2 due to a symmetry change of the ground state 
of the reference system (see Fig. usually occurs away 
from stationarity and is thus irrelevant physically. It is 
conceivable, however, that the stationary point moves to 
a point of discontinuity as a function a parameter of the 
original model. In this case the approach would generate 
an artifact which is a reminiscence of the finitcness of H' . 

So far we discussed the case U = 4 = W only. As 
a function of U the half-filled paramagnetic Hubbard 
model at T = is expected to undergo a transition from 
a metal to a Mott-Hubbard insulator £ This is marked by 
a divergence of the effective mass or, equivalently, by a 
vanishing quasi-particle weight z = 1/(1— dS(w)/cL'| tc , = o) 
as U approaches a critical value U c from below The 
result for z(U) as obtained by the use of the two-site 
reference system ( "dynamical impurity approximation" , 
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FIG. 5: U dependence of the quasi-particle weight z = 
1/(1 — c2E(w = 0)/duj) within different approximations. BR: 
Brinkman-Rice (Gutzwiller) approachi— 2-site DMFT: a non- 
variational two-site approach^- DIA, n s = 2: self-energy- 
functional approach with a reference system H' consisting of 
decoupled two-site impurity models. DIA, n s = 4: dynamical- 
impurity approximation with n s = 4 sites. DIA, n s = oo: 
DMFT limit (circles: NRGri line: ED using 8 sites). 



DIA with n s = 2) is shown in Fig. |SJ As there are less 
degrees of freedom contained in H', the approximation 
should be considered to be inferior as compared to the 
full DMFT the results of which arc in Fig. too. It is 
remarkable that the simple n s = 2-DIA (which requires 
the diagonalization of a dimer model only) yields z(U) in 
an almost quantitative agreement with the full DMFT. 

The results of the n s = 2-DIA may also be compared 
with those of the recently developed "linearized" or "two- 
site" DMFT 26 ! 28 where the Hubbard model is mapped 
onto the two-site SIAM I|1U|) by means a strongly simpli- 
fied self-consistency condition. As compared to the two- 
site DMFT, the present self-energy-functional approach 
not only represents a clear conceptual improvement but 
also improves the actual results for z{U) and for U c (see 
Fig. [SJand note that t/ c = 6 within the linearized DMFT, 
U c = 5.8450 within the DIA for n B = 2, and U c = 5.84 
and U c = 5.88 from numerical evaluations^^ of the full 
DMFT). 

The self-energy- functional approach with reference sys- 
tem of Fig. lc and small n s actually represents a 
new variant of the DMFT-exact-diagonalization method 
(ED)iiSiii As compared to previous formulations of the 
ED, the convergence with respect to n s appears to be 
faster: Compare the results for n s = 2, n s = 4 and 
n s = oo (full DMFT) in Fig. with those of Ref. [H 
More important, however, there is no need for a fit pro- 
cedure in the present approach; any arbitrariness in the 
method to find the SIAM parameters is avoided com- 
pletely. Furthermore, consistent results will be obtained 
for any finite n s while in the usual ED this can be ex- 
pected in the DMFT limit n s ^ oo only. 



V. CONCLUSIONS AND OUTLOOK 

Concluding, the proposed self-energy-functional 
method is a systematic scheme for the construction of 
new non-perturbative and consistent approximations 
for extended systems of interacting fcrmions. For 
Hubbard-type lattice models with on-site interaction 
several relations to and generalizations of existing 
approaches are obtained immediately. The numerical 
results obtained by considering a rather simple reference 
system clearly demonstrate the practicability of the 
theory. Its generality promises that the approach may 
successfully be applied also in different contexts: 

For a Hubbard-type system including M > 1 orbitals 
per site, a consistent DMFT can only be defined when 
using M baths. There is no such necessity within the 
self-energy-functional approach. While clearly the opti- 
mal local approximation requires M baths, any M' < M 
will nevertheless lead to a fully consistent approximation. 
This represents an interesting option for numerical stud- 
ies of multi-band systems. Non-local trial self-energies 
can be constructed by grouping the sites into identical 
clusters of finite size N s , switching off the inter-cluster 
hopping and treating the intra-clustcr hopping as varia- 
tional parameters. Each of the N s sites in a cluster can 
be coupled to n s — 1 additional bath sites. The relation of 
such an approach to cluster extensions of the DMFT22i2i 
and to the cluster perturbation theory22iSii2i will be the 
subject of a forthcoming paper. Extensions and appli- 
cations of the method to continuous models (inhomogc- 
neous electron gas) and to Bose systems deserve further 
investigations. 
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APPENDIX A: THE FUNCTIONAL G[E] 

For the definition of G[S], invcrtibility of the func- 
tional E[G] is required. The local invcrtibility of 
S[G] is controlled by the Jacobian T a pi a i p(iu) , iio') = 
6'E a p(i<jj)/6G a 'p'(iu>'). The two-particle self-energjsi T = 
5H/SG can be assumed to be non-singular in general. 

For a further analysis we need the following lemma: 
Consider the interaction U to be fixed. Then two dif- 
ferent Green's functions Gi and G2 must result from 
two different sets of one-particle parameters and t' 2 , 
respectively. The proof is straightforward: Consider 
the high-frequency expansion of the Green's function 
G a0 {u) = J2n=i M af) w ~ n - The coefficients are given by 
the moments = J duo w"(-l/vr) Im G a p{uj + t0+) = 

([C n c a , cp] + ) with CO = [0,H]— Using the symmetry 



7 



U a (3 7 s = Uf3 a s 7 one has: 
G a p{uj) = 5 a B — 

LO 

+ 0(lo~ 3 ) (Al) 

Now 7^ t' 2 implies the ui~ 2 coefBcients to be different 
because 

1 f°° 1 
(cUp) = -- duj lmG f3a {u + iQ+). (A2) 

Consequently, we must have Gi ^ G2. 

The lemma shows that the relation t' <-> G is one- 
to-one. Consequently, we can write S[G] = S[t'[G]] 
and T = S'S/SG = <JE/<5V • 5t'/SG with a non-singular 
Jacobian St'/SG. Hence, a singular T = S'S/SG implies 
a singular 81 /St' and vice versa. However, SI /St' is just 
the "projector" in the Euler equation Q. We conclude 
that local non-invertibility of the functional £ [G] at G = 
G(t') is indicated by 81[t'}/dt' n = with t' n = t' • n for a 
certain "direction" n in the space of hopping paramtcrs. 
For such a direction, the Euler equation J7J) would be 
satisfied trivially. 

Referring to the present numerical results, one can 
state that generally the projector 511 /St' is non-singular 
in fact, as has been expected. There is one exception, 
however, namely points in the hopping-parameter space 
where one or more bath sites are decoupled from the 
rest of the system (vanishing hybridization V) . Here the 
one-particle energy of a decoupled bath site can be var- 
ied without changing the trial self-energy. Even for this 
case, however, there are no formal difficulties with the in- 
verse functional G [£] : To ensure the local invertibility of 
£ [G] , one simply has to restrict the space of variational 



parameters t' by excluding the one-particle energies of 
the decoupled bath sites, i.e. one has to focus on the 
physically relevant parameters. This implies a respective 
restriction of the space of t'-representable Green's func- 
tions G(t') and ensures the local invertibility of £[G] 
on the restricted domain. Similarly, a restriction of the 
t' space becomes necessary to ensure the local invert- 
ibility of £[G] in case of a system where the self-energy 
is trivial (as e.g. for a model of spinless fcrmions with 
nearest-neighbor Coulomb interaction in the limit of infi- 
nite spatial dimensions where the self-energy is given by 
the Hartree term, cf. Ref . I35I) . 

Finally, it should be mentioned that generally the func- 
tional £[G] cannot be inverted globally. Consider, for 
example, the Hubbard model on the infinite-dimensional 
hypercubic lattice with nearest-neighbor hopping t at 
half-filling. Due to manifest particle-hole symmetry, a 
sign change of the hopping 1 1— s- — t leaves the (local) self- 
energy invariant but transforms (the non-local elements 
of) the Green's function G. We conclude that G[£] in 
general cannot be defined uniquely. Due to this non- 
uniqueness, and also due to non-linearity, there may be 
more than a single solution of Eq. ©. However, this does 
not cause any problem since for any £ satisfying J^J wc 
have: 

G[£] = (Go 1 - S)- 1 
=> £ = S(G ( 7 1 -S)- 1 

=> £ = £[G] and G = (G 1 - S)" 1 .(A3) 

This means that £ is given by the (formal) sum of all 
skeleton diagrams built up by a propagator G which is 
constructed via the Dyson equation from the same £ in 
turn. Hence, any stationary point should be regarded as 
a physically meaningful solution. Among different phys- 
ical solutions (corresponding e.g. to different phases) the 
minimum grand potential selects the stable one. 
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